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We present a theory of magnetization dynamics driven by spin-polarized current in terms of the 
quantum master equation. In the spin coherent state representation, the master equation becomes 
a Fokker-Planck equation, which naturally includes the spin transfer and quantum fluctuation. 
The current electron scattering state is correlated to the magnet quantum states, giving rise to 
quantum correction to the electron transport properties in the usual semiclassical theory. In the 
large spin limit, the magnetization dynamics is shown to obey the Hamilton- Jacobi equation or the 
Hamiltonian canonical equations. 
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I. INTRODUCTION 



The theory of open quantum systems, which has been 
greatly developed in the past several decades, plays a 
critical role in the understanding and control of the dy- 
namics of quantum systems that are coupled to the sur- 
rounding environments.— ~— The basic issues such as dissi- 
pation, decoherence, measurement, and noise source, etc. 
of the systems are usually investigated in the open quan- 
tum system framework. This theory has wide applica- 
tions in the fields including quantum optics^— ultracold 
atoms^ and quantum information and computation^ On 
the other hand, in micromagnetics& 7 - and spintronics^ 
the magnetization dynamics is commonly treated as clas- 
sical even though the control and dissipation parame- 
ters are couched as from quantum sources. While these 
methods have been highly successful in simulating mag- 
netization reversal and spin-torque driven magnetization 
dynamics^ there is the question whether there are quan- 
tum effects not exhibited by these classical treatments, 
when the magnet is in the mesoscopic range, of 10 3 — 10 7 
spins, between the molecular magnets and the macro- 
scopic magnets as defined by Ref. [ll|. Addressing this 
question may not only pave the way for the future tech- 
nology developments, but also broaden our vision and 
deepen our understanding of the emerging mesoscopic 
quantum world between the well established microscopic 
and macroscopic ones. In this paper, we present a the- 
ory of a single domain mesoscopic magnet as a member of 
the family of open quantum systems for its spin-current- 
driven dynamics. 

When the spin-polarized current passes through the 
ferromagnetism (FM) layer, the spin angular momentum 
of the current electrons is transferred to the FM layer 
and thus rotate the magnetizatio n 12 ' 13 . This so-called 
"spin transfer torque (STT)" effect has now become the 
most important method to control the magnetization dy- 
namics in the nano-scale structures, where the conven- 
tional Oersted field generated by the electric current is 
less practical^. Numerous magnetoelectronics devices 
have been proposed and fabricated based on the STT- 
driven magnetization dynamics^^ In spite of its great 



success and importance, the fundamental physics of STT 
in the standard theory is semiclassical. The magneti- 
zation dynamics is described by the classical Landau- 
Lifshitz-Gilbert (LLC) equation^, while the STT terms 
in the LLG equation is obtained from the quantum scat- 
tering of the spin current electrons by the classical poten- 
tial of the magnetizatio n 12 ' 13 . This semiclassical picture 
is expected to break down as the magnet is further minia- 
turized to the mesoscopic regime. Furthermore, the STT 
has been used to stimulate and control the spin waves^, 
or magnons. Therefore, a more sophisticated investiga- 
tion of the STT in the full quantum picture is necessary 
in order to adapt the new developments in the field. In 
a previous study, we shown that the continuous scatter- 
ings between the quantum macrospin state of a magnet 
and spin-polarized electrons in a simple model simulation 
not only induce the STT effect but also generate quan- 
tum fluctuations due to the quantum disentanglement 
process*^ In this paper, we will show that the quantum 
macrospin scattering model we exploited before is exactly 
solvable, and will investigate the magnetization dynamics 
from the full quantum picture by applying the standard 
theoretical techniques for open quantum systems to this 
model. 



II. QUANTUM MACROSPIN MODEL 

In parallel with the original study of STT in the semi- 
classical pictur o 12 ' 13 , we consider the motion of a single- 
domain magnet driven by the spin-polarized current. 
However, the magnet here is not described by the clas- 
sical magnetization vector M, but is represented by the 
quantum operators of the total spin angular momentum 
J. The spin-polarized electrons are injected along the 
x-direction in sequence independent of one another, and 
interact with the magnet located at x = through the 
exchange interaction. The model Hamiltonian for each 
electron interacting with the magnet is^ 
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where the first term is the kinetic term of a single current 
electron, and the second term is the interaction between 
the electrons and the magnet; Jo and J are the unit and 
total spin operators of the magnet, and s is the electron 
spin operator; the parameters Ao and A are the spin- 
independent and spin-dependent interaction strength re- 
spectively, which are used in the semiclassical model if 
the operators Jo and J are replaced by their mean field 
approximation according to the correspondence princi- 
ple. Note that the Hamiltonian is that of a free magnet 
without the external magnetic field and anisotropic crys- 
tal field. This treatment will simplify the calculations 
below without invalidating of the general conclusions. 

The STT effect originates from the elementary 
cntanglc-discntangle processes of the scattering states 
between the magnet and the spin-polarized electrons*^ 
These scattering states are deduced from the scattering 
matrix S of the Hamiltonian ((T|). Unlike the scatter- 
ing matrix in the semiclassical picture, which is defined 
only in the Hilbert space of the electron, the scattering 
matrix S in this full quantum picture is defined on the 
larger Hilbert space including both the magnet and the 
electron (see Appendix A) , which gives the STT directly 
and more informations compared with the semiclassical 
case. 



III. QUANTUM DYNAMICS EQUATIONS FOR 
MAGNET 

A. Quantum Master Equation 

Consider the scattering of the magnet spins by an in- 
jected electron as uncorrelated. The incoming compos- 
ite state of the magnet and the current electron is as- 
sumed to be a product of their respective density ma- 
trices pf n and pf n . After scattering, the outgoing states 
of the whole system, p out = Sp^ n <S> pf n S', as a result of 
the unitary scattering matrix 5, has a degree of entan- 
glement. Next, the surrounding environment decoheres 
the entangled state into a joint probability distribution 
of the possible magnet states and the corresponding elec- 
tron states. Properties of the resultant magnet state or 
the electron state are characterized by their respective 
reduced density p^ ut or electron p^ ut from tracing over 
the degrees of freedom of the other component in p ut- 
The correlated quantum dynamics of the magnet and the 
electron injected in sequence in the spin-polarized cur- 
rent drives the time evolution of the magnet state p J and 
the magnetization-dependent electron transport proper- 
ties in the electron density matrix p e . While the mean 
magnetization dynamics is within reach of the semiclas- 
sical picture, the magnetization fluctuation is given only 
by the full quantum treatment followed here without ad- 
ditional stochastic assumptions. 

For a theory of dynamics of the ferromagnet as an open 
system, we treat the current electrons as the equivalent 
of the environment. The Kraus operator— of the magnet 



is defined in terms of the scattering matrix S as the evo- 
lution operator of each encounter with a current electron, 

lCk,s;k>, s > = (k,s\S\k', s'), (2) 

with a specific basis set {\k, s}} of an incoming electron 
state of wave vector k and spin up or down state s = ±. 
We have adopted the simple model (TTJ) for the dynamics 
of the rigid macro-spin states {| J, m}} with the total spin 
number J and the z component quantum number m and 
leave the effects of spin waves for future study. Then, 
the current electron kinetic energy is conserved and the 
Kraus operators are non-zero only if k and k' are on the 
same energy shell, given by 

fCk,s;±k,s = (£ ± + S(l z , )Ck,-s;±k,s = C J s , (3) 

where the coefficients £ and £ are functions of the basic 
parameters Ao,A and fc, J. The first operator ICk. s -.±k.s 
with the same spin index s comes from scattering without 
spin transfer, and the second operator lCk,- S ;±k.s with a 
change in s represents spin transfer. These Kraus oper- 
ators are functions of the macro-spin Jo and J (see Ap- 
pendix B). In the semiclassical picture, these Kraus op- 
erators will reduce to scalars representing effective fields 
for the magnetization dynamics. 

In a scattering event, let the initial state of the cur- 
rent electron be given by the density matrix pf n = 
J2s s' fs.s'(k)\k, s)(k, s'\. This simple form may be ex- 
tended to account for a wave vector distribution or quan- 
tum coherence between different wave vectors, but will 
not be exploited here to keep the exposition simple. With 
the above Kraus operators, the quantum map of the mag- 
net state from pf n to p^ ut in the scattering is 

Pout= J2 fs,s'(k)£±k,s'';k,sp( n (1C±k,s»;k,s'?. (4) 

If the spin-polarized current is considered as a sequence 
of electrons injected at equal time interval r (a measure 
of the inverse current), the continuous evolution of the 
magnetic density matrix driven by Eq. Q , with a coarse 
graining of a time scale much larger than r, yields the 
quantum dynamics of the magnet governed by the master 
equation, 

|/(t) = Vo(i) + S(f).T(i)], (5) 

where S = Tr[<xp? n ] is the Bloch vector of the spin- 
polarized current electrons, u being the Pauli matrices. 
The operators 7o and T are polynomial functions of Jo 
and J (sec Appendix C). 7o corresponds to the unpolar- 
ized part of the current which causes fluctuations of the 
magnet motion without a net spin transfer effect. On 
the other hand, T is induced by the electron spin polar- 
ization, giving rise to both STT and the magnetization 
fluctuations. 

Note that the master equation (O is an exact solu- 
tion from the model ([T|) for arbitrary J. Thus, Eq. §5§ 
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can be applied to molecular magnets of small J. This is 
in contrast with the simulation of the quantum stochas- 
tic dynamics of the magnet in a previous studyi^ which 
used the same scattering matrix but required the large 
J condition to keep the approximation of the magnetic 
quantum state as a spin coherent state after scattering. 



B. Fokker-Planck Equation 

To facilitate computation in the large J case and, more 
importantly, to study the quantum-classical crossover of 
the magnetization dynamics, we put the master equa- 
tion |0) in the spin coherent state ^-representation^ 
analogous to the boson case. The chosen basis set is the 
overcomplete and non-orthogonal states {|J, f2)}, where 
| J, 0) is the spin coherent state of total spin J in the 
direction of f2 = (9,$). The density matrix p J in this 
representation is p J (t) = J dilPj(il, t)\ J, £2)(J, f2|, and 
the spin operators J take the form of the differential 
operators.— Substituting these expressions of p J and J 
into Eq. (|5|) with some algebraic manipulations, we ob- 
tain the Fokker-Planck equation for V j (see Appendix 
D), 



-V- (TVj) + V 2 {VPj), 



(6) 



where the unit vector m points in the direction of the 
macrospin Q = (6, <]>), the drift vector T = A(in x S) x 
m+Smx S contains the two well-known terms of STT,— & 
the diffusion coefficient T> = A(l — m • S)/ (2 J + 1) orig- 
inates from the quantum fluctuation generated by the 
scattering.— The parameters A and B are functions of the 
parameters £ and £ in the Kraus operators ([3]), namely, 
A = (2J + 1)\(\ 2 /t, B = 2%[£*C]/t (3 denoting the 
imaginary part of) which can be determined from the 
basic parameters of the Hamiltonian (TTJ. 

The quasi-probability distribution function Vj in 
Eq. © is different from the one considered in the classi- 
cal theory. Its value could be negative in some situations 
because Vj describes the quantum state of the magnet as 
faithfully as the density matrix p J . As shown in Eq. ([5]), 
the STT terms naturally arise from the open quantum 
dynamics of the magnet in the presence of the contin- 
uous scatterings by the spin-polarized electrons. Unlike 
the semiclassical picture, where the STT terms are indi- 
rectly obtained from the current electron spin polariza- 
tion after potential scattering, in the quantum case, the 
STT terms follows directly from the full quantum scatter- 
ing. Furthermore, the full quantum treatment also gives 
the quantum fluctuations accompanying the spin trans- 
fer processes, which does not exist in the semiclassical 
treatment. 

The diffusion coefficient V expression shows depen- 
dence on the relative angle between the magnet and the 
electron spin, with its maximal (minimal) value when m 
and S are anti-parallel (parallel) . This coincides with our 
previous simulation^ that the quantum magnetization 



fluctuation is first enhanced and then suppressed during 
the STT-driven magnet switching, which is observed in 
a recent experiment^ For the special case where the in- 
jected electrons are fully unpolarized (S = 0), T> is still 
a non-zero constant which again suggests that the unpo- 
larized current can cause quantum magnetization fluctu- 
ations without net spin transfer. The steady solution of 
Eq. ^ is a constant, which means a uniform distribution 
function in the spin coherent state space and the magne- 
tization will vanish on average. By contrast, the semiclas- 
sical theory predicts only the zero spin torque but no dif- 
fusion dynamics for the magnet. This STT-induced mag- 
netization fluctuation becomes dominant over the ther- 
mal magnetization fluctuation at low temperatures. The 
crossover temperature is estimated by comparing the dif- 
fusion coefficient V in Eq. ([5]) with the one for thermal 
magnetization fluctuation , 19 ' 20 
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where M is the magnetic moment of the magnet, a g 
the Gilbert damping coefficient, 7 9 the gyromagnetic ra- 
tio, ks the Boltzmann constant, and T the temperature. 
Since |M| = lg Jh and |C| 2 ~ 0(1/ J 2 ) for J > 1, wc 
obtain a g KsT ~ H/Jt, in agreement with the quantum 
noise estimated 

The Vj distribution as the solution of the Fokker- 
Planck equation © gives the exact quantum dynamics 
of the magnet based on the quantum macrospin model 
(UJ). The expectation values of any observable physics 
qualities, such as the magnetization and its fluctuations, 
can be calculated from the Vj distribution. An exam- 
ple is demonstrated in subsection D. In the semiclassical 
picture, the STT is defined on the mean field level of 
the magnetization, and quantum correlation between the 
magnetization states does not exist. The quantum the- 
ory includes the quantum correlation and is applicable to 
any possible exotic quantum states of the magnet in the 
mesoscopic or microscopic regime. 



C. WKB Approximation for Large J 

In the large J regime, we demonstrate that the solution 
of Eq. © leads to classical behavior. The expressions for 
A and B suggest that T - 0(1/ J) and V ~ 0(1/ J 2 ), 
thus the diffusion term is smaller than the drift term in 
Eq. (JHJ) by the order of magnitude 0(1/ J). Then, the 
WKB approximation is applied for large Substitut- 
ing Pj(m,t) = e ~JW(m,t) into Eq ^ p and k ce pi n g tnc 

terms to the leading order of 1/J, leads to the Hamilton- 
Jacobi equation for W(m,t), 



dW 



T VW + JV(VWf 



0. 



(8) 



Thus, the STT-driven magnet obeys the canonical dy- 
namics in the classical limit, and the function W plays 
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the role of action. For the constant spin-polarized current 
case, the corresponding Hamiltonian canonical equations 
in the spherical polar coordinates are 



d<3 
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dP$' 


dt 


~ <9$' 



(9) 



Here, the Hamiltonian function is explicitly written as 

% = J e P Q + J^P 9 + JVPl + ^-Pl. (10) 
sin 6 sin 9 

with the generalized momentum components, Pq = OqW 
and P$ = d$W , and the STT components Te and T$. 
The equations for O and <E> in ((9]) show that two more 
terms, which originate from the diffusion term in Eq. ([5]), 
contribute to the magnetic dynamics in addition to the 
STT terms even in the classical limit. Eq. ([9]) and (p~0|) 
give a more complete description of the classical magne- 
tization dynamics than the semiclassical STT theory. 




FIG. 1. (color online). Distribution function Vj for the nano- 
magnet at four different time t. (a) t = 0; (b) t — 0.3 tjv; (c) 
t = 0.5 tN\ (d) t — tjv- The simulation parameters are set as: 
J = 10 4 , A+ = 1.3 V, A_ = 0.1 V, d = 3 nm, N = 1.5 X 10 5 , 
k = 13.6 nm~\ S = (0, 0, 1). A 200x200 lattice in 9-$ plane 
is used in the simulations. 



D. Numerical Example 

Here we demonstrate how to apply the approach de- 
veloped above to the STT-driven quantum dynamics of 
a magnet. We consider a magnet with J = 10 4 , and the 
initial distribution function obeys the Boltzmann distri- 
bution, i.e. Vj(rh, 0) = Ce~ E l kBT , where the energy of 
the magnet in a magnetic field B is = — M • B, C is 
the normalization factor, and the magnetic moment is 
M = -fghJm. We set the temperature as T = 1 K, and 
the magnitude of the magnetic field B = 0.05 T. The di- 
rection of B is chosen that maximum value of Vj locates 
at the angle fio = (2.8, 1.0). The initial distribution is 
show in the Fig. [Ha). Then we apply a spin current pulse, 
which includes 1.5 x 10 5 electrons. The Bloch vector of 
the electron spin is S = (0,0,1), and the wavevector is 
k = 13.6 nm -1 . To calculate the scattering matrix, the 
parameters Ao and A in Eq. ([1} are estimated for a mag- 
net with effective potentials A + = 1.3 V, A_ = 0.1 V 
and thickness d = 3 nm. 

In the simulations, we have used the method of 
characteristics^ to solve the Hamilton- Jacobi equation 
(|SJ), which gives the time-evolution of W(m,t) and then 
the distribution function Vj(m,t) = e - JW ( m ^\ Direct 
solution of the Fokkcr-Planck equation or the mas- 
ter equation ([S]) is also practical but not explored here. 
The time is measured in tjq — Nr. The distribution 
functions Vj at t = 0.3tN,0.5t^,tpj are shown in Fig. Q] 
(b)(c)(d) respectively. In order to keep the normaliza- 
tion ofVj, it is renormalized after every 500 scatterings. 
Note the different scales for 9 in Fig. Q] We found that 
Vj is first expanded and then compressed, and its center 
moves from (2.8,1.0) to (0.05,4.2), which show the effect 
of spin current on the distribution function. 
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FIG. 2. (color online). Time-evolution of J M=a; ,j / ,z and its 
fluctuation 5in= x ,y,z (x: solid line; y: dash line; z: dotted 
line.) caused by the spin transfer torque obtained from the 
probability distribution function (thick line) and the quantum 
trajectory method (thin line). 



The mean value of the macrospin J and its fluctuations 
are calculated from Vj as 



J^=»,»,«(*)= / dnv{n,t){j,n\^\j,n), 



sr m , z (t) = / dnr(n,t)(j,n\6^\j,n). 



The results are shown in Fig. [21 where the mean magne- 
tization is switched by the STT, and the magnetization 
fluctuations are first enhanced and finally suppressed. 
Comparison with the results obtained from the quantum 
trajectory method in Ref.[l5| gives reasonable agreement. 
The approach developed here together with the quantum 
trajectory method^ fits in the toolbox for micromagnet- 
ics simulations with the added value of accounting for 
relevant quantum effects. 
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IV. ELECTRON DENSITY MATRIX AFTER 
SCATTERING 



Finally, we briefly discuss the electron states after the 
scattering, which contain the informations about the elec- 
tric current and current noise, etc. After tracing over 
the degrees of freedom for the magnet in the total den- 
sity matrix p out , the reduced density matrix of the elec- 
tron is p e out = Trj[Sp e ln <g> p{ n S^]. This expression is 
a generalization of the electron potential scattering in 
the semiclassical picture 2 ^. Here, the transformation of 
the electron state is no longer unitary due to the recoil 
of the magnet. For instance, we assume that one elec- 
tron with wavevector k > and spin-polarized vector 
S = (0,0,1) is injected, i.e., p\ n = \k,+)(k,+\, and take 
the 'P-representation for the quantum states of the mag- 
net. With the scattering matrix S, we obtain 



Pout 



= J2 \k',s')(k",s"\ 

k' ,s' ,k" ,s" 

x / dnp(n)<n|(x: fc « ja « ;fc , + )t/c fc i 1 ./. fcl+ |n>. (11) 



For the model ([T]), the terms in Eq. (| 1 1 1) are non-zero 
only if k' = ±fc and k" = ±fc. The electron scattering 
state is correlated with the quantum state of the magnet. 
The scattering formalism in the semiclassical picture will 
be reproduced if the Kraus operators are replaced by the 
corresponding scattering matrix elements there. As the 
magnets are miniaturized further— and the quantum de- 
scription becomes necessary, the semiclassical scattering 
formalism for the electron transport will break down, and 
Eq. (fTTj) and its generalizations should be exploited as the 
new starting point. 



V. CONCLUSION 

In conclusion, the STT-driven magnetization dynam- 
ics has been investigated by treating the magnet as an 
open quantum system in the exactly solvable quantum 
macrospin model. A set of dynamical equations is es- 
tablished and the quantum-classical connection is made. 
The full quantum picture here provides a unified and 
complete description of the magnetization dynamics and 
electron transport, and further explorations of the quan- 
tum physics in spintronics along this line is expected. 
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Appendix A: Scattering Matrix 

First, we calculate the scattering matrix S of the model 
Hamiltonian 



H = ~d 2 x +5{x) (A J + As -J 



(Al) 



Considering that one electron in state \ij)f n ) is injected 
along the x-direction and the initial quantum state of 
the magnet is \4>i n }, then the incoming state |*m) of the 
whole system before scattering is the product state of 
|VfJ and \iPi n ), i.e., |* in ) = \^ n ) ® After scat- 

tering, the outgoing state \SS> out ) will be |* ou t) = S\$ in ). 
The scattering matrix S are determined by the boundary 
conditions at x = 0, 

(*m + *o«t)U=0- = (*m + *o«t)U=0+, 



I H(^ in + y out )dx = e (<£ m + * 
Jo- Jq- 



out)dx{A2) 



where e is the total energy of the whole system. 

The scattering problem (|A2[) is simplified by utilizing 
the symmetries of the model Hamiltonian {AT}. First, 
the kinetic energy of the electron is conserved during the 
scattering process. Thus the scattering matrix elements 
of S is non-zero only if the absolute values of the incoming 
and outgoing wavevectors of the electron are the same. 
Second, the operators (s+J) 2 and? 2 -|-J 2 are commutative 
with the Hamiltonian (|A1|) . and their eigenstates \J,p) 



arc given as 



(s + J) 2 \J,p) 
(s z +3 z )\J,p) 



J(J+l)\J,p), 



with J = J ± i and p = — J,..., J. Choosing the basis 
set {\k\J, p)}, the scattering problem (|A2j) reduced to 
a set of 5-potential scattering equations, and gives the 
scattering matrix S in this representation^. Then after 
a representation transformation with the help of Clebsch- 
Gorden coefficients, we obtain the scattering matrix S in 
the basis set {|fc, s; J, to}}, which takes a block form 



/ ■ 



S 





o s kttl 



V o o '• 

and the form of each block Sk tfi is 



(A3) 



k,fi 



( 



\ 



k.LL 

k.u k.u 



k,U 



k,u 



r k,u 
t 



^k,u 



k.u k.u k.u 



(A4) 



k,u k.u k.u k,u 



We acknowledge the support of this work by the U. 
S. Army Research Office under contract number ARO- 
MURI W911NF-08-2-0032 and NSF ECCS-1202583. 



Here, the element t s k s u (r^* ) is the transmission (re- 
flection) probability amplitude from \k,s';J,p — ^s'} to 
\k, s; J, p — (| — k, s; J, p — 5 s ))- The spin transfer 
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is related to those elements with s ^ s'. The explicit 
expressions for the matrix elements are 



J.+ + —IT? 71 2 I — 171 7 _ -2 

tkjl = cos VJ,+ e ' + cos + cos?7./,_e sin 

+ + ■/ • —177 7 _j_ 2 i — Z71 7 _ ■ 2 \ 

r fc,M = -H sln lJ,+ e cos + sin ?7 j, _e M " sm aj lfl ), 
h,u = cos77j, + e~ ,,w,+ sin 2 aj i(1 + cos r)j,-e~ lrlJ -- cos 2 aj i(1 , 



i(sin??j i+ e " ?J '+ sin 2 qj :M + sinT^-e 



-ill j — 2 



COS CtJ.i 



h^t = (cos rjj >+ e' lvj - + — cos77j i _e _I ' ,J ' _ ) sina,/ iM cosqj iM , 
rj~+ = -i(sin77j i+ e~ !,,J '+ - sin rjj i -e~ lvj ^ ) sin a JiM cos aj i(1 , 
f fe,M = (cos77j i+ e _l,?J - + — cos77j,_e _1,,J ^ ) sina,/ iM cosqj iM , 
= -i(sin?7j i+ e~ ! ' ,J ' + — sin r]j t -e~ lvj ^ ) sin aj i(1 cos 

Here, the phase shifts arc given as rjj,± = tan -1 A i ± , 
with the effective potentials Aj,+ = ^( JAo + JA),Aj,_ = 
i [JAo — (J + 1) A]. The Clebsch-Gordan coefficients 

cosa,/ M and sinaj !M are given as coaaj ifl = J '^j+^ 'i 



2,7+1 



Appendix B: Kraus Operators 

Now we express the Kraus operators ICk.s-.k'.s' in the 
basis set {|J, m)} based on the scattering matrix S ob- 
tained above. The block form of S means that ICk,s;k',s' is 
non-zero only if k and k' have the same absolute values. 
For example, we have 

- (k,+\S\k,+) 



V 











tl 



■ t 



(Bl) 



which is a (2J + l)-dimcnsion diagonal matrix. With the 
Clebsch-Gordan coefficients, the diagonal elements are 
rewritten as 



where 



c = 



j+i 



2,7+1 



k,r 



cos r?j, + e 



""7J,+ 



+ 



J+l 



2,7+1 



sm?7j i+ e 



-"Jj,+ 



2J + 1 
J 



cos rjj t -e 



+ 



27 + 1 



sm ?7j,-e 



- 1 »7J,- 



+ 



2,7+1 
1 



2J + 1 



(cosi7j,+e %VJ <+ — cos r]j t -e iVj > ) 
(sin r)j t J r e~ lT,J - + — sin»7j,-e * VJ '~). 



Considering the matrix form of the angular momentum 
operator J z in the basis set {|J, m}}, the matrix (|B1[) 



shows that the Kraus operator Kk,+;k,+ is just 
/C fe ,+ ;fe , + = (£+i)X + cX, 

where Jo is the unit matrix. 

Similarly, the other Kraus operators are written in the 
compact form as 

lCk,s;±k,s = (£ ± 7:)Jo + S C-^z, fck,-s;±k,s = CJs- 



Appendix C: Quantum Master Equation 

The master equation ([5]) in the main text is obtained by 
substituting the Kraus operators ([3]) into Eq. (0| there. 
The calculations are straightforward, and yield the ex- 
plicit expressions for the operators 7o and T as 

% = (\e - \) P j + ici 2 (X/x + + h.c, 

% = 2^cVX + ici 2 (X/j+ - J+/X) + h.c, 

T y = 2tt*p J % + z|C| 2 (J+p J X -Xp J X) + h.c, 

% = 2^cVX + ICI 2 (J+/X -X P J J+) + h.c. 

The terms containing a single angular momentum oper- 
ators in T can be interpreted as a field-like torque, while 
the terms including two angular momentum operators 
the Slonczewski-type torque and the quantum fluctua- 
tions. This becomes clearer in the spin coherent state 
representation. 



Appendix D: The Fokker-Planck Equation 

Here we explain the derivation of the Fokker-Planck 
equation (6) in the paper. In the spin coherent state rep- 
resentation {|J, il)}, the density matrix pj is expressed 
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(Dl) 



dnvj{n)\j,n)(j,n\ 



Using S = (a, (3, 7), and substituting the expression (|D1[) 
into the master equation (5) in the paper, and utilizing 
the differential forms of the operators^ J^| J, ft) (J, ^|Jj 
= 0, +, — , z), we derive the differential equation for 
Vj(ty as 



dVj 
dt 



1 d(-T @ Vj 



sine 
1 



sin 6 99 



'smB 



1 d(-T*Pj) 
sin 6 9$ 

I d 2 {VVj) 



80 



sin 2 9 



(D2) 



Here, 



Te = A(a cos 8 cos $ + (3 cos O sin $ — 7 sin 9) 
+ B(a sin <£> - cos <f>) , 



7 



T$ = A(— a sin $ + /3 cos $) 

+ B(a cos O cos $ + /3 cos 6 sin $ — 7sin0), 
A 

V = (1 — a sin 9 cos $ — 8 sin 6 sin $ — 7 cos 6), 

2J+1 V P ' ; 

with the coefficients A = (2J + l)^f-,B = 2a( f C) . Fur- 
ther analysis shows that Te and T$ are the components 
of the spin transfer torque T = „4(m x S) x m + Bin x S 
in the spherical coordinates, where the unit vector m 



denotes the direction of the macrospin. Replacement of 
the differential operators in spherical coordinates by the 
divergence operator V and Laplace operator V 2 reduces 
Eq. (|D2[) to the simple form of the Fokkcr-Planck equa- 
tion (FPE) 

^Pj(nM) = -V ■(TPj) + V 2 (VVj). (D3) 



* Present address: Department of Physics, University of 

Hong Kong 
t lsham@ucsd.edu 

1 H. J. Carmichael, Statistical methods in quantum optics 1: 
Master equations and Fokker-Planck equations (Springer, 
Berlin, 1999). 

2 H. J. Carmichael, Statistical methods in quantum optics 2: 
Nonclassical fields (Springer, Berlin, 2008). 

3 H.-P. Breuer and F. Petruccione, The Theory of Open 
Quantum Systems (Oxford University Press, 2002). 

4 F. T. Arecchi, E. Courtens, R. Gilmore, and H. Thomas, 
Phys. Rev. A 6, 2211 (1972). 

5 C. Weedbrook, S. Pirandola, R. Garca-Patron, N. J. Cerf, 
T. C. Ralph, J. H. Shapiro, and S. Lloyd, Rev. Mod. Phys. 
84, 621 (2012). 

6 A. Aharoni, Introduction to the Theory of Ferromagnetism 
(Clarendon Press, Oxford, 1996). 

7 H. Suhl, Relaxation processes in micromagnetics (Oxford 
University Press, 2007). 

8 D. Ralph and M. Stiles, Journal of Magnetism and Mag- 
netic Materials 320, 1190 (2008). 

9 A. Brataas, A. D. Kent, and H. Ohno, Nat. Mater 11, 372 
(2012). 

J. Z. Sun, Phys. Rev. B 62, 570 (2000). 

1 W. Wernsdorfer, Advances in Chemical Physics 118, 99 
(2007). 

2 J. C. Slonczewski, Journal of Magnetism and Magnetic Ma- 
terials 159, LI (1996). 



L. Berger, Phys. Rev. B 54, 9353 (1996). 

M. Madami, S. Bonetti, G. Consolo, S. Tacchi, G. Carlotti, 

G. Gubbiotti, F. B. Mancoff, M. A. Yar, and J. Akerman, 
Nature Nanotech. 6, 635 (2011). 

Y. Wang and L. J. Sham, Phys. Rev. B 85, 092403 (2012). 

K. Krause, States, Effects, And Operations Fundamental 

Notions Of Quantum Theory, edited by A. Bohm, J. Dol- 

lard, and W. Wootters (Springer- Verlag, Berlin, 1983). 

L. M. Narducci, C. M. Bowden, V. Bluemel, G. P. Gar- 

razana, and R. A. Tuft, Phys. Rev. A 11, 973 (1975). 

V. E. Demidov, S. Urazhdin, E. R. J. Edwards, M. D. 

Stiles, R. D. McMichael, and S. O. Demokritov, Phys. Rev. 

Lett. 107, 107204 (2011). 

W. F. Brown, Phys. Rev. 130, 1677 (1963). 

Z. Li and S. Zhang, Phys. Rev. B 69, 134416 (2004). 

H. Risken, The Fokker-Planck Equation (Springer- Verlag, 
Berlin, 1984). 

R. Courant and D. Hilbert, Methods of Mathematical 
Physics (Interscience, New York, 1953). 
J. Foros, A. Brataas, Y. Tserkovnyak, and G. E.W. Bauer, 
Phys. Rev. Lett. 95, 016601 (2005). 

J. K. W. Yang, Y. Chen, T. Huang, H. Duan, N. Thiya- 
garajah, H. K. Hui, S. H. Leong, and V. Ng, Nanotechnol- 
ogy 22, 385301 (2011). 



